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^ ABSTRACT 

We investigate threshold pion photoproduction in the frame'work of heavy baryon chiral 
/\ ' perturbation theory. We give the expansion of the electric dipole amplitude -E0+ to three 

c^ I orders in n, the ratio of the pion to nucleon mass, and show that it is slowly converging. 

We argue that this observable is not a good testing ground for the chiral dynamics of 
QCD. In contrast, we exhibit new and fastly converging low-energy theorems in the 
P-waves which should be used to constrain the data analysis. We also discuss the 
importance of polarization observables to accurately pin down certain multipoles and 
give predictions for the reaction 771 -^ n^n. 
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I. INTRODUCTION 

Over the last few years, much interest has been focused on pion photo- and elec- 
troproduction off nucleons. In particular, new accurate data for the processes 7p -^ ir^p 
and 7*p — > ir^p close to production threshold have become available [1,2]. These have led 
to many theoretical investigations. In particular, in refs.[3] baryon chiral perturbation 
theory was used to give a model-independent description of the pertinent differential 
cross sections, amplitudes and so on. In these papers, the nucleons were treated as 
fully relativistic fields which leads to complications in the power counting underlying 
the effective field theory [4]. These can be overcome by a clever choice of the spin-1/2 
fields being velocity-dependent which allows to transform the baryon mass term into 
a string of 1/m (here, m denotes the nucleon mass) suppressed interactions [5,6]. In 
this formulation, called heavy baryon chiral perturbation theory (HBCHPT), there is 
a strict correspondence between the loop expansion and the expansion in small mo- 
menta and quark masses. Our motivation to come back to the topic of threshold pion 
photoproduction off protons and neutrons is fivefold. First, the calculation in the rela- 
tivistic framework indicated that the expansion of the electric dipole amplitude i?o+ is 
indeed slowly converging as a function of /U = Mt^/iji, with M^^ denoting the pion mass. 
However, in that approach the complete expression of order ^^ could not be given. In 
HBCHPT, if one calculates within the one-loop approximation but to next-to-next- 
to-leading order g^ [q denotes any small momentum or mass), the first three terms in 
the chiral expansion of -^0+ can be given. Second, in the relativistic formulation [3] we 
found that the differential cross sections were not well described due to an essentially 
energy-independent i?o+ of — 1.3- 10~^/M^+. Also, in the P-waves the relative strength 
of the Ml- multipole came out too large as compared to the Mi+. Here, we will show 
that two of the P-waves are severely constrained by novel low-energy theorems and that 
the third one is completely dominated by the A(1232) resonance. These LETs were im- 
plicitely contained in the relativistic calculation but not made explicit. Furthermore, 
we give here a much better estimation of the pertinent low-energy constants based on 
the idea of resonance saturation. This in turn leads to a satisfactory description of the 
existing data in the threshold region. Third, it is obvious that the information gained 
from the total and differential cross sections is not sufficient to pin down all multipoles 
uniquely. For doing that, one has to consider polarization observables and we will dis- 
cuss some of these here. Last, but not least, we also show detailed predictions for the 
reaction 771 -^ n^n which is experimentally very difficult to access but in fact has to 
be studied for various reasons, one of them being the test of the isospin decomposition 
based on first order electromagnetism which is usually assumed in the construction of 
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the invariant amplitudes. Finally, new data for 7p -^ n^p in the threshold region and 
above taken at MAMI (Mainz) and SAL (Saskatoon) are presently being analyzed. 

The manuscript is organized as follows. In section 2, we fix our notation and define 
the pertinent observables to be discussed later. In section 3, we briefiy discuss the 
effective pion-nucleon Lagrangian underlying our calculation. Sections 4 and 5 contain 
the main results of this paper. We present the order q^ calculation for the S-wave 
including some isospin-breaking from the pion mass difference and the corresponding 
0{q^) calculation for the P-waves.^-* We then discuss the low-energy theorems (LETs) 
in the S- and P-waves followed by the presentation and discussion of the numerical 
results. There are two low-energy constants entering the expression for Eqj^{uj) (with 
u) the pion cms energy). These can either be fixed by a best fit to the available data 
or estimated by resonance exchange. Already here we would like to stress that the first 
method leads to unnaturally large numbers for these coefficients which are a refiection 
of the importance of higher loops not yet calculated. Consequently, even to this order 
there remains some appreciable theoretical uncertainty leading to the conclusion that 
the electric dipole amplitude is in fact not the best testing ground of the chiral dynamics 
of QCD as was long believed. We also present a two-parameter model for Eqj^{u) which 
simulates most of the physics in the threshold region. Furthermore, the imaginary part 
of i?o+ and its relation to the Watson final state theorem is discussed showing again 
that to this order the description of the S-wave is not too accurate. However, we also 
stress that the presently available determinations of -Eo+ close to threshold hinge on a 
few empirical points because the P-waves quickly dominate the cross section. Therefore, 
it is imperative to study polarization observables. These allow for a clean separation of 
certain multipoles and a much more accurate empirical determination of small multipoles 
like e.g. -E'i+. The summary and outlook is given in section 6 and some lengthy formulae 
are collected in the appendix. 

II. THRESHOLD PION PHOTOPRODUCTION: FORMAL ASPECTS 

In this section, we will give the formalism necessary to treat pion photoproduction 
in the threshold region. We will only be concerned with the kinematics close to threshold 
and the corresponding multipoles. We also summarize the formulae for the differential 
and the total cross sections as well as for some polarization observables. 

Consider the process 7(/c) + N{pi) -^ Tv^{q) + N{p2), with A^ denoting the nucleon 
(proton or neutron), 7 a real (/c^ = 0) photon and tt^ the neutral pion. The polarization 



^^ Due to the fast convergence in these multipoles, a more accurate calculation does not 
seem necessary. 



vector of the photon is denoted by e^. In the threshold region, the three-momentum 
q of the pion in the TiN centre-of-mass (cm) frame is smaU and vanishes at threshold. 
It is therefore advantageous to perform a multipole decomposition since at threshold 
only the S-wave survives and close to threshold one can confine oneself to S- and 
P-waves. The corresponding multipoles are called {E, M)i±, where E, M stands for 
electric and magnetic, / = 0,1,2,... the pion orbital angular momentum and the ± 
refers to the total angular momentum of the pion-nucleon system, j = I ±1/2. These 
multipoles parametrize the structure of the nucleon as probed with low energy photons. 
Consequently, the T-matrix T ■ e depends on four multipoles and takes the following 
form in the cm system 

in /v /v ^ 

T ■ e = ia ■ e{Eo+ + k ■ qPi) + ia ■ ke- qP2 + {q x k) ■ ePs (2.1) 



The quantities -Pi, 2, 3 represent the following combinations of the three P-waves, 

Ei+, Mi+ and Mi_, 

Pi = 3Ei+ + Mi+ -Ml- 

P2 = 3Ei+ - Mi+ + Mi_ (2.2) 

P3 = 2Mi+ + Mi_ 

These four amplitudes are calculable within CHPT. As will become clear later, the 
particular choice (2.2) of the P-waves is best suited for the chiral expansion and the 
physics related to it. 

We now discuss briefly the kinematics for 7p -^ n^p. The pion energy in the cm 
system is given by 

s-ml + M^o ^ E +Ml,/2m^ 

2v^ ^1 + 2E^/mp 

with s the cms energy squared and E^ = [s — rn^) / {2mp) the photon energy in the lab 
frame. At threshold, Sthr = {mp + MT^o)'^ and cuthr = "^o = Mt^o. The second threshold is 
related to the opening of the 7r+n channel at Uc = 140.11 MeV (since Sc = {'mn+M^+)'^). 
In what follows, we will take into account the pion mass difference M^± — M^o = 4.6 
MeV since it subsumes the most important isospin-breaking effects. This is discussed 
in some detail in ref.[7]. However, we do not differentiate between the proton and the 
neutron mass in the loops. Consequently, we will use M^± = ujc = 140.11 MeV to 
account for the proper location of the second threshold. The tiny error induced by this 
procedure is well within the theoretical uncertainty of our approach. 
The differential cross section can be written as 

O!^ = A + B cose + C cos" e (2.4) 

\q I ail an 



in the approximation that only S- and P-waves contribute. 9 is the cms scattering 
angle, \k\ = {s — rn^) / {2^) and \q\ = JuS^ — M^^. The energy-dependent coefficients 
A, B and C are related to the multipoles via 

I U+l -r ^1 21 -r ^1 .i\ 

B = 2Re{Eo+P*) (2.5) 

1 1 

r^ I P |2 I P |2 I P |2 

I ii 2' ^' 2' ' 

These are real for ujq < uj < uJc and complex above ujc (we do not consider here the 
tiny phase related to the direct iv^p scattering process [8] because it only shows up at 
two-loop accuracy). The total cross section follows as 

|^(7tot = 47rfA+ Jc) . (2.6) 

From the nearly forward-backward symmetric angular distributions exhibited by the 
Mainz data, one can immediately conclude that |A|, |C| ^ \B\ which means that a 
very accurate knowledge of the P-waves is mandatory to reliably extract the electric 
dipole amplitude. Therefore, different assumptions on the P-waves can lead to a rather 
different energy variation of the electric dipole amplitude i?o+(cu) in the threshold region 
[9]. Also, in most analysis it is assumed that i?i_|_ = 0. To get a handle on such small 
multipoles and to allow for a clean separation of the various real and imaginary parts, one 
has to investigate polarization observables. We will consider here the polarized photon 
asymmetry S(^), the polarized target asymmetry T{9) and the recoil polarization P{9). 
These are given by 

E(^)=rsin^(|P3|'-|i'2n 

T{e) = 2rim((Eo+ + cos^Pi)(P3 - ^^2)*) (2.7) 

P(0) = 2rim((Eo+ + cos^Pi)(P2 + P3)*) 



with 



\q\sme [ d(T\ 

: • (2-7a) 
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In fact, P{0) and T{9) allow for a direct determination of Im i?o+ above cVc since the P- 
waves are essentially real in the threshold region (as we will show later). This concludes 
the necessary formalism. 



III. EFFECTIVE LAGRANGIAN 

In this section, we will briefly discuss the chiral effective Lagrangian underlying our 
calculation. To explore in a systematic fashion the consequences of spontaneous and 
explicit chiral symmetry breaking of QCD, we make use of baryon chiral perturbation 
theory (in the heavy mass formulation) (HBCHPT). The nucleon mass is considered 
large compared to typical momenta in the system. This allows to decompose the nucleon 
Dirac spinor into "large" (H) and "small" (h) components "^{x) = e~'^^'"'^{H{x)+h{x)} 
with v^ the nucleon four- velocity, v'^ = 1, and the velocity eigenfields are defined via 
i/H = H and i/h = —h. Eliminating the "small" component field h (which generates 
1/m corrections), the leading order chiral tcN Lagrangian reads 

6^i^H{ivD + gAS-u)H (3.1) 

Here the pions are collected in a SU(2) matrix-valued field U{x) 

U{x) = 4 [A/F2-7f(x)2 + if ■ 7f (x)] (3.2) 

r 

with F the pion decay constant in the chiral limit and the so-called cr-model gauge 
has been chosen which is of particular convenience for our calculations in the nucleon 
sector. In eq.(3.1) D^ = d^ + V^ denotes the nucleon chiral covariant derivative, S^ is 
a covariant generalization of the Pauli spin vector, qa — 1-26 the nucleon axial vector 
coupling constant (formally the one in the chiral limit) and u^ = iu^V ^Uu^ ^ with 
u = VU and V^ the covariant derivative acting on the pion fields. To leading order, 
0{q) one has to calculate tree diagrams from 

4A^ + ^Tr{V'^t/V^t/t+x+} , x±=u^XU^±ux^u (3.3) 

where the second term is the lowest order mesonic chiral effective Lagrangian, the non- 
linear a-model coupled to external sources. The scalar source x is proportional to the 
quark mass matrix Ai. Beyond leading order, the effective Lagrangian takes the form 

r„- £(1) + £(2) , £(3) , £(4) , £(2,4) .3 4N 

C^^j*^' contain 1/m corrections and counterterms. The a priori unknown coefficients of 
these counterterms are the so-called low energy constants. For the calculation to order 
g^, one has to consider tree diagrams with insertions from >C^^' as well as one-loop 



diagrams with insertions from C^^ . The terms from C^j^ of relevance to the problem 
at hand are 



£[5^5^]((l + ^„)/i + ^^Tr/+)}if 



(3.5) 



where f+^u = e{u^Qu + uQu^)F^^j, Q = (1 + T3)/2 and F^i, is the canonical photon 
field strength tensor. Note that k,s,Vj the isoscalar and isovector anomalous magnetic 
moments of the nucleon in the chiral limit, are combinations of the low-energy constants 
cq and C7 discussed in [7] . We note here that the other low-energy constants of order 
g^, which are called 01,2,3,4, do not contribute at all to our final results for the chiral 
expansion of the S- and P-waves (to the order we are working). There are two terms 
from /^jriv which enter [10], 



4'ir = hUi^ vxe^^^'P Tr (/+^.t/p) + 622 {ivxe^^'"' S, /|^, vD + h.c.) jif + . . . (3. 



6) 



with f^fj^i, = f+nu ~ (1/2) Tr/_|_^j^. Although the coefficients 614 and 622 are infinite [10], 
in the case at hand only the finite combination 6p ~ 4 614 + g^^ 622 contributes and the 
low-energy constant bp enters the P-wave P3 in the following way: 

Pi\uj)=ebpuj\q\ . (3.7) 

Note that for the other multipoles i?o+ and Pi, 2 one can not construct any gauge and 
Lorentz invariant chirally symmetric counter term at 0{q^) respecting the discrete sym- 
metries P, C, and T. Finally, the minimal set of terms from C^j^ which give a contribution 
to Eo^{uj) (for the proton) is of the form 

£Sv = ^{ -47rF a2 S^ v"" /+ {v Dv u) - 2mF ai S^ v" /+ X- W + . . . (3.8) 

and show up in Eo+ in the following way, 

E^\{u) = e ai{X)u M^ + e a2iX) u^ . (3.9) 

These counter terms are necessary to absorb the divergences generated by the loops at 
order g^. Here, A is the scale introduced via the dimensional regularization. Ideally, 
this scale dependence is cancelled from the one in the corresponding loop contribution. 
If one, however, estimates the contact terms via resonance saturation, one is left with 
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a small but spurious scale dependence as detailed in ref.[13]. We also note that the 
low-energy constants bp, ai and a2 are a priori different for the proton and the neutron. 
The dots in eqs. (3.6,3.8) stand for the respective chiral vertices accompanied by the 
low-energy constants. We omit the superscript '(P?'^)' here since from the context it is 
obvious which reaction we are considering. There are furthermore many terms in C^^.^ 
which come from the l/m-expansion of the relativistic ttN Dirac Lagrangian. These 
generate the 1/m-expanded Born contributions to the S- and P-waves given in the 
appendix. Finally, we note that we will work in the Coulomb gauge cq = since in that 
gauge one has not direct lowest-order photon-nucleon coupling (which is proportional 
to e ■ v) and thus many diagrams vanish. 

IV. QCD ANALYSIS 

In this section, we briefly outline the 0{q'^) calculations for the P-waves and the 
0{q'^) for i?o+- The explicit expressions are relegated to the appendix. Next, we estimate 
the low-energy constants which enter the expressions for Eo+ and P3. We then turn 
to the discussion of the low-energy theorem (LET) for i?o+ and exhibit novel LETs for 
the P-waves. This constitutes one of the major results of this investigation. We end 
this section with a short discussion of the imaginary part of -Eo+ and its relation to the 
Fermi- Watson theorem. 



IV.l. CHIRAL EXPANSION TO ORDER 0{q 



3^ 



We wish to calculate the T-matrix elements to order q . Since the photon polar- 
ization vector counts as 0{q), one gets the S- and P-wave multipoles with an accuracy 
of 0{q'^). In the Coulomb gauge, we have tree contributions from C\^^ to Eq^ and Pi 
{i = 1, 2, 3) of the type 

where /2, /s and Qi are dimensionless functions of their arguments. In the one-loop 
diagrams, we include the pion mass difference to account for the most important isospin- 
breaking effect. £'o+(cu) is given by the triangle and rescattering diagrams (see fig.la,b) 
and the other two diagrams shown in that flgure contribute to Pi^2{^)- To this order, 
there is no loop contribution to ^3(0;). The pion mass difference in the loops leads to 
the cusp in Eqj^ of the square-root type 




Eo+('^) = coefficient ■ Wl 5- + ... (4.2) 



where the eUipsis stands for polynomial pieces (in uj). We will come back to the coeffi- 
cient multiplying the square root when we discuss the final state theorem. At this order, 
there is no mass and coupling constant renormalization due to the vanishing of the tree 
couplings. Finally, there are also finite counter terms at this order. The anomalous 
magnetic moment of the proton (neutron) enters the expressions for i?o+ and -Pi, 2- In 
addition, there is the finite coefficient bp defined in eq.(3.7) contributing to P3. The 
estimate of its numerical value will be given in section 4.3. 

IV.2. CHIRAL EXPANSION TO ORDER 0{q^) 

For the electric dipole amplitude, we consider one more order. This is motivated by 
the fact that the relativistic calculation of ref. [3] points towards the importance of higher 
orders. On more general grounds, we remind the reader of the well-known fact that in 
S-wave observables it is often mandatory to go beyond 0{q'^). However, we should also 
stress already at this point that a calculation up-to-and-including 0{q'^) might not 
be sufficiently accurate. Furthermore, we do not differentiate here between the neutral 
and charged pion masses and use only M^+ to get the proper cusp effect at this order. 
In the polynomial pieces (of the loops) , this infiicts a theoretical error proportional to 
{uj — M^+)/M^+ ~ 4% which is smaller than the uncertainty from the determination of 
the low-energy constants. The dominant unitarity (cusp) effect is already accounted for 
at order g^, compare eq.(4.2). At order g^, there are tree diagrams which lead to terms 
of the type uj^ jra^ and there are many loop graphs. The latter can be categorized as 
follows: (i) loop graphs with one vertex from Cy^]^ and one vertex from C\l^^ (ii) with 
vertices from L\1^ but a nucleon propagator from L\1^ and (iii) the 1/m corrections to 
the loops calculated at order (^ . Of course, many of the loop diagrams account for mass 
and coupling constant renormalization, 

(m, g^^, Kp^n, F, M) -^ (m, ^f^Ar, Kp,^, F^, M-,,) . (4.3) 

Finally, there are the novel counterterms which contribute to i?o+ as given in eq.(3.9). 
Therefore, the expressions for E'o+ and -Pi, 2, 3 calculated to order q^ and g^, respectively, 
take the generic form 



E,^{u) = E^ri^) + E^Y'^^^ioo) + E,\-'--^{u;) + E^,X{ 



UJ 



P,{uj) = Pf'^'^iuj) + P,^ -'°°P(c^) , z = 1, 2 (4.4) 



where 'Born' subsumes the nucleon-pole and anomalous magnetic moment contribu- 
tions.^^ The exphcit expressions for the various terms appearing in eq.(4.4) can be 
found in the appendix. Here, we just remark that the Born and counterterm contribu- 
tions are real and that the loop contributions are complex for u > Uc- We also point 
out again that the leading terms for i?o+ and Pi^2 appear at the same chiral power, 
namely 0{q^). Indeed, the P-waves are proportional to \q\, but not to \q\ \k\ as usually 
assumed. However, since \k\ only varies by four per cent from o^thr to Uc, this makes 
no visible effect. Note furthermore that P3 is essentially given by the contact term P3* 
since the Born contribution to this multipole is very small. It is also instructive to notice 
that the magnetic part of the nucleon Born terms is dominated by the Mi_|_ and Mi_ 
multipoles with 2Mi_|_ + Mi_ ~ whereas static A exchange leads approximatively to 
Pi = — P2 — -P3/4, i.e. 2Mi_|_ -f- Ml- is much larger than Mi+ — Mi_. This is, of course, 
a particular feature of the threshold region, further up in energy the Mi_|_ multipole 
quickly becomes dominant. 

IV.3. ESTIMATION OF THE LOW-ENERGY CONSTANTS 

The most difficult task is to pin down the values of the low-energy constants (LECs) 
ai, 02 and hp. We concentrate on the reaction 7p -^ n^p and mention the necessary 
modifications for the case 771 -^ n^n in the end of this section. We will follow two 
approaches here. In the first one, we will use the available total and differential cross 
section data to fix these coefficients. However, as we will see, this leads to unnaturally 
large values of the constants ai and 02 because they subsume the effects of higher loops 
not yet calculated. Therefore and secondly, we will estimate the LECs making use of 
the resonance saturation principle [11]. This can be formulated as follows. Consider me- 
son resonances (M=V,A,S,P) and baryonic excitations (A^* = A(1232), A^*(1440), . . .) 
chirally coupled to the Goldstone bosons (collected in U) and the matter fields (A^). 
Integrating out the meson and nucleon excitations, 

f [dM] [dN*] exp i I dx £eff [U, M, A^, A^*] = e^p i f dx C^s [U, N] (4.5) 

one is left with a string of higher dimensional operators contributing to £eff [U, N] in 
a manifestly chirally invariant manner and with coefficients given entirely in terms of 
resonance masses and coupling constants of these resonance fields to the Goldstone 
bosons. A specific example for the baryon sector is discussed in ref.[12]. For the case at 



^) This decomposition facilitates the comparison with the existing literature but is not 
a consequence of the chiral expansion. 
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hand, we have mesonic and baryonic contributions. As aheady discussed in [3], there is 
t-channel vector meson exchange, here the p(770) and the c<;(782). These contribute to 
ai,2 as follows {V 



p'^ + io) 



aX 



ASTT^mFS 



4.45 GeV" 



tti 



2 V 



-1.78 GeV" 



(4.6) 



where we have used QpN = fi^/S, gu;N = 3(7/2, Mp = M^ = y^gFj^, g = 5.8545, Kp = 6, 
K^ = and G^p^ = g/{167r'^F^) = (1/3)0^^;^ together with m = 938.28 MeV and 
Ftt = 93 MeV. Notice that the values of these coupling constants are somewhat simplified 
and given in part by the gauged Wess-Zumino action. However, as it will turn out, the 
vector meson contribution to the P-wave P3 is rather small and in the S-wave, there 
are uncertainties due to higher order effects so that the accuracy given by these values 
is sufficient. We could as well take the vector dominance value of k^^j = —0.12 (and so 
on) without any noticeable change in the numerical results to be discussed later. The 
sign of the vector meson contribution is fixed from the sign of the triangle anomaly. The 
P-wave contribution is found as 



.V 



647r3F3 



3.13 GeV"^ 



(4.7) 



The result given in eq.(4.7) is not affected by the tensor coupling in agreement with the 
considerations presented in ref.[3]. 

From the baryon sector, the by far largest contribution comes from the A(1232) 
resonance. It is mandatory to consider the A as a fully relativistic spin-3/2 field before 
integrating it out. Therefore, the contribution to the various low-energy constants will 
depend on the two AA^7 coupling constants gi and (72, the irNA coupling (/ttata = 
(3/\/2)(7^7V and the off-shell parameters X, Y and Z (see also the discussion in ref.[13]). 
This procedure leads to 



tti 



+ 



Cgi 
6 

Cg2 



ttn 



+ 



Cgi 
12 

Cg^ 
16 



—2m\ + 2m^m + m^ 
niA — m 

m{2X + 1) +mA 

lOm^ — TmAfn — Sm^ 
niA — m 

m{2X + 1)(1 - 6Z) + 2mA (4XZ + X + Z) 



+ 2m{2Y -Z- 2YZ) - 2mA{Y + Z + AY Z) 



(4.8a) 



+ 2m(5y + 2Z- 2YZ) + 10mA(l" + Z + AY Z) 



(4.86) 
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with C = Qt^n / {Q\/2Tim^m\) = 0.40 GeV~^ and Qt^n = 13.4. Throughout, we use 
the Goldberger-Treiman relation to fix g^, i.e. g^ = girNFT^/m. In what foUows, we 
will adopt two strategies. First, we keep (71 = (72 = 5 fixed and vary X, Y, Z in the 
ranges given in ref. [14] and then also allow to vary gi and (72 within the ranges given in 
[14]. As it will turn out, the results are fairly insensitive to the variation of gi and g2- 
Inspection of eqs.(4.8a,b) reveals a rather large uncertainty infiicted from the relatively 
poor knowledge of the off-shell parameters. We will come back to this when we discuss 
the numerical results. Furthermore, the A contributes prominently to 6p, 



,A Cgim 



2ml + ruAm - m^ ^ ^ ^ ^ ^^ +2m^iY + Z + AY Z) 

TTlA — m 



(4.9) 



Notice that to the order we are working, 6p receives no contribution proportional to g2- 
To get an idea about the size of 6p, we use the static isobar model (which involves no 
off-shell parameters) and find 



^A.static _ 9l97TN{mA - m) 



^A,static ^ ^ ^ yiy.rA.in.zX-n., ^ ^^.S GeV^ (4.10) 



which is considerably larger than the vector meson contribution, eq.(4.8). For the neu- 
tron, the only difference is the sign of the p^-meson contribution. We get 

ar'" = lar'^(z = l,2); b^'^ = U^'^ . (4.11) 

In this case, we have no data to fit these LECs and must use the resonance exchange 
estimates. 

IV.4. LOW-ENERGY THEOREMS (LETs) FOR THE S- AND P-WAVES 

We consider here the expansion in fx = Mt^/tti of -Eo+i -Pi and P2 at threshold. 
Since P3 is completely dominated by the contact term proportional to 6p, a similar 
expansion for this combination of the P-wave multipoles does not make sense. First, 
however, let us briefiy state what is meant by a LET following ref. [15]. We consider as 
a LET of order q^ a general prediction of CHPT to 0{q'^). General prediction means a 
strict consequence of the Standard Model depending on some low-energy constants like 
Ftt, m^gA, Kp, ■ ■ •, but without any model assumption for these parameters. This gives 
a precise prescription for obtaining higher-order corrections to the leading order LETs 
which can e.g. be obtained from current algebra. 
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First, we study the electric dipole amplitude. For that, we work in the isospin limit 
niu = iTid and to first order in the electromagnetic coupling constant. The //-expansion 
takes the form 



-^0+,thr 



Stttti 



^{Cl+ ^iCl+ ^i^Cl +o(^^)} 
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(4.12a) 

At present, the LETs given by (4.12) do not have too much predictive power since the 
only way to determine the LECs ai(A) and a2(A) are the threshold data from Mainz [1] 
(for the proton). Using the best fit values (see section 5), one finds a\{m) +a^{m) = 2.7 
GeV""^ leading to^) 



-^0+,thr 



-3.45 ( 1 - 1.26 + 0.55 +...)■ 10~^/M^, 



-1.0 ■ 10"V^7r+ (4.13) 



where the ellipsis stands for terms of order (U^ and higher. Setting ai{m) + 03 (m) = 0, 
the 0.55 would read 0.64 and the corresponding value for Eq_^ ^.^^ = —1.33 ■ 10~^/M7r+. 
For the neutron, we find (using the same values for the LECs) E^^ ^^^ = 3.64 ( 1 — 
0.29 + ...)■ 10~^/M^+ = 2.59 • 10~^/M^+ which shows a better convergence since the 
term of order ^ is absent and thus the contribution from the triangle diagram appears 
already at lowest order. We note that the electric dipole amplitude for yr*^ production 
off neutrons is sizeable and of opposite sign to the one for production off protons. The 
lesson to be learned is that the fu. expansion of i?o+ converges very slowly (as already 



^' Notice that the individual values of ai{m) and 03 (m) are much larger than their sum 
so that this determination is afflicted with a substantial uncertainty. Even the sum ai + 02 
is strongly affected if one chooses either the charged or the neutral pion mass in eq.(3.9). 
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anticipated in ref.[3]) and thus one has at least to go one order higher before one can 
make a reasonably accurate theoretical prediction. This clearly shows that the electric 
dipole amplitude is not a good testing ground for the chiral dynamics of QCD. Such a 
behaviour is, however, not too surprising. We remind the reader of similar large higher 
order S-wave effects in the scalar form factor of the nucleon [16] or the scalar form factor 
of the pion [17] (just to name two such cases). To further tighten the prediction on £'0+, 
it is conceivable that one has to perform a dispersive analysis supplemented with CHPT 
constraints to get a handle on the higher orders. It is again important to stress that the 
LETs for the electric dipole amplitude have been derived in the exact isospin limit. It 
is not meaningful to compare the number (4.13) with the data, since isospin breaking 
and other higher order effects are substantial. 

We now turn to the P-waves. From the explicit formulae given in the appendix we 
derive the LETs for the slope of Pi and P2 at threshold (for the proton) 



1 TD ^9nN J ^ . 
■pq--Pl,thr = -Z ^ <^1 + Kp + ^ 

\q I 



2 487r 
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Sum? \ 






+ 0{fi' 
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" 127r 
with Up = 1.793 and similarly for the reaction 'jn -^ n^n, 
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(4.14a) 
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(4.146) 



with Kn = —1.913 the anomalous magnetic moment of the neutron. These are examples 
of quickly converging fj, expansions,"^-* 

-^Pl^^^ = 0.512 (1 - 0.062) GeV"2 = 0.480 GeV"^ 

T^^Plthr = -0.512 (1 - 0.0008) GeV-2 - -0.512 GeV-^ 

^Pl^^j^^ = 0.351 (1 - 0.020) GeV"2 = 0.344 GeV"^ 

^^2"thr = -0.351 (1 + 0.107) GeV"^ = -0.389 GeV"^ 



(4.15) 



^^ To be precise, we mean that the leading term is much bigger than the first correction 
in contrast to what happens in the electric dipole amplitude. 
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From these expressions, a few interesting observations can be made. First, we note that 
the multipole Ei^ is not exactly zero since in that case one would have Pi = —P2- 
Commonly, this multipole is set to zero when one analyzes the threshold data. We 
will come back to this small multipole when we discuss the polarization observables. 
Second, these P-wave LETs help to constrain the existing fits to the threshold region 
[9] . They favor the solution which leads to a strong energy-dependence of the electric 
dipole amplitude. If we pull out by hand a factor |A; |^-*, these LETs translate into 

Pi%r = l0.3\k\\q\10-yM^^ P2V = -10.91^11^-110-3/^1^3+ (4.16) 

to be compared e.g. with the value of Pi = (8.8±0.6)|/c ||g llQ-^/M^.,. given by Drechsel 
and Tiator [18]. Third, we also would like to stress that the corresponding Pi and P2 
of the relativistic calculation [3] agree quite nicely with the LET (remember that in the 
relativistic formulation some higher order terms are included). For example, at Ej = 151 
MeV, the LET predicts Pi = 2.47 and P2 = -2.48 while the P-wave multipoles of ref.[3] 
lead to Pi = 2.43 and P2 = -2.60 (all in units of 10-^/M^+) (for -fp -^ n^p). These 
novel P-wave LETs should be tested more accurately and they can also serve to constrain 
the empirical analysis. It is amusing to note that this is a good testing ground for chiral 
dynamics in contrast to common folklore. 

IV.5. RELATION TO THE FERMI-WATSON THEOREM 

Here, we wish to elaborate briefly on the imaginary part of the electric dipole 
amplitude. By virtue of the Fermi- Watson theorem, it is related to the tiN scattering 
phases via (see e.g. [20]) 

-1 r\ 

liaEll^ = Re^i^J tan(5i) + -Re^J^/^^ tan(5i) + -ReE^^^^ tan(53) (4.17) 

with 5i^3 the -kN S-wave phases for total isospin 1/2 and 3/2, respectively. Close to 
threshold, we can approximate the tan(52/) by a2i \q\ (in the respective channels) and 
also drop the term proportional to a"*" Re Eqj^ which is a factor of 200 smaller than 
V2a-ReEJ+", i.e. 



li^Ell^ c^V2a- M^ReEZt'' J^-l . (4.18) 



^^ We stress again that this is not the correct threshold behaviour of the P-wave 
multipoles. 
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Therefore, the strength of the imaginary part of i?o+ in the threshold region is governed 
by the product V2a- M^ReE^^"^ = 3.7 ■ 10-^/M^+. This is in fact the coefficient 
which we did not write exphcitely in eq.(4.2). The order g^ calculation leads to 

^(l-|i)^2.T.10-VM,.. (4.19) 

which shows that the strength of Im Eo+ is underestimated by approximately 30%. As 
already noted a couple of times, this indicates that for an accurate description of the 
electric dipole amplitude even in the threshold region one has to go beyond order g^. 
Of course, we should also stress that the imaginary part is in any case less accurately 
calculated. Here, the first contributions to Re i?o+ are of order q^ whereas the corre- 
sponding imaginary part starts at 0{q^). Finally, we wish to point out why the phase 
related to the direct w^p scattering process only appears at two-loop accuracy making 
use of the Fermi- Watson theorem. Above the tt^p but below the yr+n threshold, we have 
Im Eqj^^ = Re Eq_^ ■ ta,n d{n^p). In the T-matrix, a~^ starts at order q^ and Re Eq_^ is 
0{q). Furthermore, there is an additional factor q from the relation tan^^^ ^ = a~^ \q\ so 
that the imaginary part starts out at order g^, i.e. is 0{q^) in the full amplitude which 
is a two-loop effect. 

V. RESULTS AND DISCUSSION 

In this section, we will discuss first a two-parameter model which describes most of 
the physics in the threshold region and draw some general conclusions from it. We then 
turn to the detailed numerical investigation of the reaction 'jp -^ n^p making use of the 
existing Mainz data [1]. We will perform a free fit to the three low-energy constants 
ai, a2 and bp but also one constrained by resonance exchange considerations. We argue 
that the latter one is presumably more realistic. We also show the pertinent P-waves 
and discuss polarization observables. We then turn to some predictions for neutral pion 
production off neutrons, with the LECs estimated via resonance exchange. 

V.l. A REALISTIC TWO-PARAMETER MODEL FOR Eo+iu) 

Tree diagrams and resonance exchanges lead to an almost constant Re i?o+ in the 
threshold region, 144.7 MeV < E^ < 160 MeV. The unitarity corrections due to the 
opening of the yr+n channel at 6.8 MeV above threshold have a square-root behaviour 
as discussed in section 4.1. This claim is further substantiated by the fact that in the 
isospin limit one finds essentially no energy dependence in Re Eq^ [3] but only after 
inclusion of the pion mass difference (and, to a lesser extent, the proton-neutron mass 
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difference) a strong energy dependence develops [19]. To a good approximation, we can 
therefore parametrize the electric dipole amplitude in the threshold region as 

E^+{uj) = -a-h^l-{u/u,Y (5.1) 

which immediately leads to a square-root type behaviour for Im -Eo+, 

luiEo^{uj) =h^{uj/uJcY -lQ{u-uJc) (5.2) 

which means that the coefficient h is constrained by the Fermi-Watson theorem (cf. 
section 4.5). We now fit the Mainz data with this form for i?o+ and the P-waves as 
given by the chiral expansion, i.e. we have three parameters, namely a, h and 6p, 
to fit 126 data points (total and differential cross sections). We find using standard 
minimization procedures 

a = (0.28 ± 0.07) ■ 10"^^^+, 

6= (4.62±0.49)-10-VM^+, (5-3) 

bp = (15.64 ± 0.25) GeV-^ 

Several remarks on these numbers are in order. First, the value for bp is close to the 
static A exchange estimate eq.(4.10), i.e. the multipole P3 is completely dominated by 
A exchange even close to threshold. ^-^ Second, the fitted value for b is somewhat larger 
than what one would get from the Fermi-Watson theorem, bpw = 3.7 (in canonical 
units of 10~^/M^+). The source of this discrepancy is two-fold. First, in the deriva- 
tion of Im i?o+ from the Fermi- Watson theorem we assumed exact isospin symmetry 
and made the further approximation that the phase shift is simply the product of the 
scattering length times the momentum. The result obtained was, however, applied to a 
situation involving some isospin breaking. Second, if the remeasured threshold data lead 
to somewhat smaller values of Eo+ in the threshold region, this difference of 20% would 
diminish. Furthermore, the rather simple but physically motivated form eqs.(5.1,5.2) 
leads to a good fit with a x^ /datum of 1.89. At the respective thresholds, this gives 

ReEQ+{uo) = -1.52 ■ 10"^^^+ ReEo+{uc) = -0.28 ■ 10"^^^+ (5.4) 

which translates into a difference of SEq^ = Eo^{uJc) — Eqj^{ujo) = 1.24 ■ 10~^/M^+ 
between the two thresholds. Due to the refiection properties of (5.1), this also means that 



^' The sign of the multipole P3 is determined from existing multipole analyses at some- 
what higher energies, see e.g. refs. [24] . Of course, only in the polarization observables to 
be discussed later this sign plays a role. 
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the imaginary part of Eqj^ at iVji = 2uJc — ujq = 145.25 MeV (equivalent to E^ = 158.27 
MeV) should be 

luiEo+{uR) = 1.24 ■ 10"VM^+ • (5.5) 

Stated differently, a strong variation of the real part between the yr^p and the yr+n 
thresholds reflects itself in a rapid growth of the imaginary part and vice versa. This 
stringent constraint rooted in dispersion theory has not yet been discussed in the various 
examinations of the energy-dependence of the electric dipole amplitude as given by 
the Mainz data. If one uses a square-root behaviour of the imaginary part in the 
threshold region, the value of Bergstrom [9], Im i?o+(180 MeV) = 2.0 translates into 
Im Eoj^(ujr) = 1.0 and the four values below 182 MeV given by Miillensiefen [20] lead 
to Im Eq^{ujr) = (l.OiO.l) (all in canonical units). These numbers are consistent with 
a direct calculation of Im Eo^{ujr) = 3.7 ■ ^uj\Iuj\ — 1 = 1.0, and they indicate that 
the variation of Re i?o+ between the yr^p and the yr+n thresholds is indeed less strong as 
commonly believed. We furthermore stress that the relative smallness of the parameter 
a = 0.3 indicates that there have to be large corrections to the Born result of 2.3 (in 
canonical units). Only with the inclusion of loop diagrams, here the triangle graph 
and its crossed partner, it is possible to understand such large corrections to a. This 
can be considered a success of CHPT. A last important point is the following. In the 
approximation (5.1), Re Eo^{uj) is strictly constant for uj > Uc- At present, the data are 
not accurate enough to clearly differentiate between a constant or slowly varying energy 
dependence above Uc- We have therefore added a linear term of the type c (1 —uj/uJc) to 
eq.(5.1) and redone the fitting. The values for a and h are somewhat changed leading to 
ReEo+(c^o) = -1.57-10-^^7^+ and ReEQ+{ujc) = -0A3-10-^/M^+, not very different 
form eq.(5.4) with a comparable x^/datum of 1.89. However, the one-cr uncertainties 
on a and b are considerably larger for this type of fit. This means that a reshuffling 
between the linear and the square-root term is possible (using the existing data). We 
will come back to this when we discuss the fit with the counter terms proportional to 
ai(A) and a2(A). We end this section by stressing again our believe that the pertinent 
ingredients of the threshold behaviour of the electric dipole amplitude are indeed given 
by the form eqs.(5.1,5.2). 

V.2. CHPT RESULTS FOR 7p -^ n^p 

We now turn to the discussion of the results making full use of the formalism 
outlined in section 4. First, we consider ai(A), a2(A) and bp as completely unconstrained 
parameters and use the best fit to the Mainz total and differential cross section data 
to determine their values. This will be called the "free fit" in what follows. Second, 
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we vary these coefficients witliin tlie bounds given by tfie resonance excliange picture as 
discussed in section 4.3. Tfiis means in particular tfiat we vary tfie off-sfiell parameters 
X, Y and Z (for fixed gi = g2 = 5). In principle, one sfiould also vary tfie vector meson 
couplings witfiin some bounds but as discussed before, tfie uncertainty witli respect to 
tlie A parameters is by far larger and we thus use a fixed vector meson contribution to 
the various LECs. We will, however, also present a fit in which gi and (72 are allowed 
to vary within their bounds. We already note here that the results are essentially 
indistinguishable from the fit with gi = g2 = 5. This procedure will be coined the 
"resonance fit". With M^+ ~ 140.11 MeV, the yr+n threshold is located at its physical 
value, E^ = 151.43 MeV. 

First, we show results for the free fit. We find 

ai(m) = (-55.45 ± 3.34) GeV"'^ , 

a^{m) = (58.15 ± 3.14) GeV"^ , (5.6) 

bp = (15.80 ± 0.23) GeV"^ 

Notice that the value for bp is in excellent agreement with the resonance exchange 
estimate, b^ + b^ = (12.3 + 3.1) GeV-^ = 15.4 GeV'^, using eqs.(4.7,4.10). Such 
a value is essentially a consequence of the bell-shaped differential cross sections for 
E^ > 150 MeV. The result of the free fit for the differential cross sections is shown in 
Fig. 2 by the solid line (the fit to the data at E-^ = 156.1 MeV is not exhibited) and 
similarly in Fig. 3 for the total cross section.^-* This fit has a x^/datum of 1.88. If one 
multiplies the values of ai(?n) and a2(m) by M^, one notices indeed that their individual 
contributions to i?o+ at threshold are of the order ±5.6-10~^/M^+ which is considerably 
larger than their sum. The corresponding real part of the electric dipole amplitude is 
shown in fig. 4a, with 

ReEo+iuo) = -1.56 ■ 10"^^^+ ReEo+{uc) = -0.32 ■ 10"^^^+ (5-7) 

not very different from the constant plus square-root fit discussed in the previous section. 
However, after the Tr+n threshold. Re i?o+(cu) rises in contrast to the two-parameter 
model. We believe that this is an artefact of the strong energy dependence from the 
polynomial contact terms proportional to ai^2 with their large coefficients. This is remi- 
niscent of the fit we discussed before when we added a linearly growing term to eq.(5.1). 
We interpret the unnaturally large values of ai(?7i) and a2{m) as a signal of the im- 
portance of higher loop effects not accounted for by our order g^ calculation. After all. 



^) Notice that the few Saclay data were not used in the fitting procedure. Including 
them would not alter any of our conclusions. 
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from the discussion of the unitarity corrections leading to the cusp effect it is rather 
obvious that one can not expect a strong energy dependence due to some polynomial 
terms (in u). This is exactly what happens in this free fit - the unconstrained param- 
eters try to make up for some higher order effects. Clearly, the situation would be 
much more satisfactory if one could determine the LECs ai,2 from some other reac- 
tion. The imaginary part of the electric dipole amplitude (cf. fig. 4b) is not affected 
by such uncertainties. It shows the expected square-root type rise and stays below 
0.7 ■ 10~^/M^+ for Ej < 160 MeV as already elaborated on in section 4.5. In Fig. 5, 
we exhibit the conventional P-wave mulitpoles Mi_|_, Mi_ and i?i+. These show the 
empirically expected pattern Mi_|_ > — Mi_ ^ E'l^ but still Pi ~ — P2 and the A is 
most visible in P3. We note that the P-waves are improved compared to the relativistic 
0{q^) calculation [3]. It is furthermore important to stress that the imaginary parts of 
these P-waves are tiny (they increase with \q\^) and that consequently the cusp at the 
TT+n threshold is not visible. Isospin breaking effects are also small, typically of the size 
((M^+ - M^o)/M^+)3/2 ~ 0.6%. The energy-dependence of Pi,2/\q\ in the threshold 
region ujq < uj < ujr is very weak, i.e. these reduced P-waves stay constant on the 
level of 2%. These observations are at the heart of the usefulness of the P-wave LETs, 
eqs.(4.14). 

Before discussing the polarization observables, let us consider the resonance fit. 
First, we keep (71 = (72 = 5 fixed. ^^ We find as best values with a x^/datum of 2.02 

X ^ 2.24 ± 1.87 , Y = 0.13 ± 0.52 , Z = 0.28 ± 0.75 . (5.8) 

This corresponds to af{m) = 1.3 GeV'^, ai^im) = 2.7 GeV'^ and bp = 15.9 GeV-^. 
As already stated, the magnitude of ai^2(?^) is considerably smaller as in the case of the 
free fit for the reasons discussed. The corresponding differential and total cross sections 
are shown in Figs.l and 2 as the dashed lines. They are very similar to the free fit, the 
sole exception being the first two MeV above the n^p threshold. The resonance fit leads 
to a smaller i?o+,thr and weaker energy dependence as shown in Fig. 4a. Specifically, we 
have 

ReEo+{uo) = -1.16 ■ 10"^^^+ ReEo+{uc) = -0.43 ■ 10"VM^+ (5.9) 

This means that the cusp is less pronounced. We point out, however, that the energy- 
dependence of i?o+ for the resonance fit follows closely the generic form constant plus 
square root, eq.(5.1) (with somewhat different values for a and b). The large one-cr 



Note that the empirical width r(A -^ N'j) demands gi ~ 5. 
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uncertainties on X, Y and Z signal that the presently available data are not yet accu- 
rate enough to pin down the electric dipole amplitude very tightly. In the resonance fit, 
Re Eo^{u > Uc) stays flat as does the two-parameter model discussed above. Conse- 
quently, the relation between the imaginary part at un and the difference in the real 
part between uq and Uc is fulfilled (cf fig. 4b) . The P-waves are essentially the same as 
for the free fit. If we relax the condition that gi = g2 = 5, we find a very similar fit with 
the following values: X = 0.55 ±0.32, y = 0.65 ±0.02, Z = 0.24 ±0.01, £?i =3.94 ±0.06 
and (72 = 4.49 ± 4.24. Since g2 enters only the S-wave, its value is determined within 
large uncertainties. For example, this 5-parameter fit leads to i?o+,thr = —1.17, only 
marginally different from the 3-parameter fit. 

We now turn to the polarization observables defined in eq.(2.7). In Fig. 6, we show 
E(6'), T{9) and P{0) for E^ = 153.7 MeV. Since none of these is sensitive to Re -Eo+i 
the predictions based on the free fit and on the resonance fits are essentially the same 
(compare the solid to the dashed lines in Fig. 6). From the size of the effect we conclude 
that the target asymmetry T{9) is best suited to pin down the imaginary part of the 
electric dipole amplitude. The photon asymmetry E is very sensitive to the ratio of the 
P-wave multipoles i?i+/Mi_ (for fixed Mi+).^-* However, note that our analysis gives 
a very small i?i+ in the threshold region so that the sensitivity discussed in ref.[21] is 
presumably overestimated. 

V.3. CHPT PREDICTIONS FOR 771 -^ n^n 

We now discuss predictions for neutral pion production off the neutron in the thresh- 
old region. To fix the low-energy constants, we use the resonance exchange values dis- 
cussed in section 4.3 (since no data to fit exist). Also, we set m = rrin = 939.57 MeV. 
The threshold of this reaction is at E!:^"^ = 144.66 MeV and the p7r~ channel opens at 
E^ = 148.46 MeV, i.e. at ujc = 137.86 MeV, and we choose the charged pion mass 
accordingly to account for the proper location of the second threshold. 

The resulting total cross section from threshold up to E^ = 160 MeV is shown in 
fig. 7 and four corresponding differential cross sections in fig. 8. We note that crtot rises 
quicker than in the case of the proton, this is due to the dual effects of (i) the larger (in 
magnitude) electric dipole amplitude, -E'o+" changes from 2.13 to 2.77 between threshold 
and E-f = 160 MeV (in canonical units), and (ii) the even closer proximitiy of the first 
open channel. One also notices the cusp effect. The differential cross sections are 
strongly peaked in forward direction. This can be traced back to the large and positive 
value of Re Eq^, implying a large and positive angular coefficient B (eqs.(2.4,2.5)). This 



^) This was pointed out to us by R. Beck, see also ref.[21]. 
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is very different to the case of the proton. At cUc, we find -E'o+"^ = 2.79, i.e the cusp in 
the electric dipole is of similar size as in the proton case. The polarization observables 
shown in fig. 9 (for E^ = 153.7 MeV) are strongly enhanced in backward direction which 
is due to the forward peaked da/dQcm- It is obvious that these rather distintive features 
should be tested experimentally. 

VI. SUMMARY, CONCLUSIONS AND OUTLOOK 

In this paper, we have used heavy baryon chiral perturbation theory to study the 
reactions 'jp -^ Tr^p and 'jn — * jr^n in the threshold region. This is a continuation and 
improvement on the calculations making use of relativistic baryon CHPT reported in 
refs.[3,19]. The pertinent results of this study can be summarized as follows (preliminary 
results were presented in ref.[22]): 

• In the threshold region, one can restrict oneself to the inclusion of S- and P-wave 
multipoles as defined in eq.(2.1). We have calculated the electric dipole amplitude 
to order q^ and the P-waves -Pi, 2, 3 to order q^. Besides the Born and loop contribu- 
tions, we have one finite counter term contribution to P3 and two scale-dependent 
ones to -E'o+- These counter terms can either be determined by a best fit (the 
so-called "free fit") to the threshold data (total and differential cross sections) or 
estimated from resonance exchange ("resonance fit"). In the latter case, the dom- 
inant contributions come from the A(1232) as well as the vector mesons p and iv. 
In both cases we get a good fit to the existing data. However, for the free fit, the 
two S-wave low-energy constants are unnaturally large and of opposite sign. This 
signals the importance of higher loop effects not yet accounted for. The P-wave 
low-energy constant is essentially given by A-exchange and takes a value expected 
from the static isobar model. 

• We have considered the low-energy theorem for Eqj^ and shown that the con- 
vergence in (U = M-i^/m is indeed very slow (as conjectured in ref.[3]), compare 
eq.(4.13). We conclude that this multipole is not a good testing ground for the 
chiral dynamics of QCD. In case of the neutron, the large contribution of order 
Ij? from the triangle diagram and its crossed partner appears already at leading 
order and thus the convergence for Eq^^ is much better. From our calculation, one 
expects that \Ef_^''\ > \Ef/\. 

• We have derived novel low-energy theorems for the P-waves Pi and P2 as given 
in eqs.(4.14,4.15). These are quickly converging expansions in p and they should 
be used to constrain the data analysis. In particular, the small difference in the 
magnitudes of Pi and P2 indicates a small but non-vanishing i?i+ multipole. 

21 



• We have presented a simple but realistic two-parameter model for the energy depe- 
dendence of the electric dipole amplitude in the threshold region, cf. eqs.(5.1,5.2). 
The parameter b is closely related to the strength of the imaginary part of i?o+ by 
the Fermi-Watson theorem. In fact, the energy dependence of the real part should 
reflect itself by a similar rise in the imaginary part above the yr+n threshold. This 
points towards the importance of an independent determination of the imaginary 
part at E^ = 158.3 MeV. 

• We have discussed the polarization obervables ^(6*), P{0) and T{9) and shown that 
the target asymmetry T seems to be best suited to determine Im i?o+ for the neutral 
pion production off protons. 

• Finally, we have given predictions for the total and differential cross sections as 
well as for polarization observables for the reation 'jn — ^ ir^n (with the low-energy 
constants determined from resonance exchange). These should be determined ex- 
perimentally since they serve as a further test of the chiral dynamics of QCD. 

Where do we go from here? It is imperative to improve upon the S-wave on the the- 
oretical side by either a two-loop calculation or a dispersive representation constrained 
by CHPT and on the experimental side by more accurate determinations of the total 
and differential cross section close to threshold. Only with very accurate data one is able 
to test the proposed constant plus square-root form (5.1) for the S-wave constrained 
by the Fermi- Watson theorem and the P-wave LETs. Furthermore, a new look at the 
corrections to the low-energy theorems for charged pion photoproduction seems to be 
required from the new data on tt" production [23]. We hope to come back to these 
topics in the future. 
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APPENDIX: EXPRESSIONS FOR S- AND P-WAVE MULTIPOLES 

Here, we will give explicit analytical expressions for the S- and P-wave multipoles 
i^o+j -Pi, -P25 -P3 of neutral pion photoproduction from protons and neutrons. The for- 
mulae are given in the isospin limit using only one pion mass M^r. In order to account 
for the branch point and unitarity cusp above the physical threshold, the value of Mj^ 
has to be chosen appropriately, M^+ = ujc, as explained in the text. The Born terms 
for the proton read: 

QoOTrm"* y w^ u J 

(Al) 



+ 



^ ^ ^ 87rm2 \ 5m 5mc^ ^V 2my ' ^ ^ 



Stttti^ I^ 10m 5ma; \ 2m/ 

•^ ^ ' 167rm3 ^ ^ 

The Born terms for the neutron are: 



8nm'^ 2m 






P^''''^{u) = (A.7) 

The loop contributions for the proton and the neutron differ only by a few numerical 
coefficients: 
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+ VM2 - c^2 LmI + ^[(11 + 4t^)cj^ + (1 + 2t^)M^] arcsin ^ j 



ujMf. ( 27r + arcsin — -— 



UJ 



arcsin 



M^ 



+ 



9\ 



A 3 3 a; w^ 






r2^M^ 



U) 



Mt ?>Ml 



+ 37rM^(^ - 4c^) arcsin -— + 3M^(^ + 



PrP(a;) 



647r2F3 1 3^2 



UJ arcsin 



UJ 

Ml 



(A9) 



,loop,, ,, _ e^^^m j J_^^3 _ (m2 - c^2)3/2] ^ ^^ _ y^ 



2 - M^ . ^ 

^ — uj^ arcsin — 



CO 



P"°-(-) = iiJf {3^2 [A4= - (A4^ - -^)='1 - Af. 



jloop 



P^°°P(C^) = 



(A.IO) 
(All) 



(A.12) 



The analytic continuation of these expressions above the branch point uj = Mt^ is ob- 
tained through the following substitutions: 



VM2-c^2 



-i\luj^-Ml^ arcsin > - + z In ^ (A.13) 



'tt ' 



The counterterm contributions for the protons and neutrons are: 



ESV(^) = e aTW ^^ + e <'"(A) c^ M^ 

P^\uj) = P^\uj) = 
P^\uj) = eb''p''uj\q\ 



(A14) 

(A15) 
(A.16) 
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FIGURE CAPTIONS 

Fig.l Feynman diagrams contributing at 0{q'^). The triangle (a) and rescattering (b) 

diagrams give Eoj^{uj) while (c) and (d) contribute to the P-wave multipoles Pi 

and P2- Crossed graphs are not shown. 
Fig. 2 Differential cross sections for 'jp -^ rr^p. The solid lines refer to the free fit and 

the dashed ones to the resonance fit as explained in the text. The data are from 

Mainz [1]. 

Fig. 3 Total cross section for 'jp -^ n^p. For notations, see fig. 2. The data are from Mainz 
(diamonds) and Saclay (squares) [1]. 

Fig. 4 The electric dipole amplitude for •yp -^ ir^p. (a) The real part. For notations, see 
fig. 2. In addition, the Icr-bands for the free fit are indicated by the dotted lines. 
(b) Imaginary part. 

Fig. 5 The P-wave multipoles Mi_|_, Mi_ and i?i+ for the free fit. 

Fig. 6 The polarization observables S(^), T(9) and P{0) for E^ = 153.7 MeV. For nota- 
tions, see fig. 2. 

Fig. 7 Chiral prediction for crtot(7?^ -^ n^n). The low-energy constants are estimated from 
resonance exchange. 

Fig. 8 Chiral prediction for the differential cross section for 'jn -^ n^n at E^ = 149.1 
(dashed), 151.4 (dotted), 153.7 (sohd) and 156.1 (dash-dotted) MeV. 

Fig. 9 Chiral prediction for the polarization observables E, P, T for 'yn — ^ n^n at E^ — 
153.7 MeV. 
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